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Abstract 

The general Minimal Supersymmetric Standard Model introduces 
new sources for Lepton Flavor Violation (LFV) as well as CP-violation. 
In this paper, we show that when both sources are present, the elec- 
tric dipole moment of the electron, dg, receives a contribution from the 
phase of the trilinear A-teim of staus, (j)A^- For (j)A^ = '/r/2, the value 
of de, depending on the ratios of the LFV mass elements, can range 
between zero and three orders of magnitude above the present bound. 
We show that the present bound on rules out a large portion of 
the CP-violating and the LFV parameter space which is consistent 
with the bounds on the LFV rare decays. We show that studying the 
correlation between d^ and the P-odd asymmetry in r — > 67 helps us 
to derive a more conclusive bound on (j)A^ ■ We also discuss the possi- 
bility of cancelation among the contributions of different CP-violating 
phases to de- 



1 Introduction 

As is well-known, elementary particles can possess electric dipole moments 
only if CP is violated. In the framework of the Standard Model (SM) of 
the elementary particles, the 3x3 quark mixing matrix (the CKM matrix) 
can accommodate a CP-violating phase. In fact, CP-symmetry in the meson 
system has been observed to be violated in accord with the SM. However, 
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the effect of tlie CP-violating pliase of tlie CKM matrix on tlie Electric 
Dipole Moment (EDM) of the electron, rfg, is shown to be very small [1] and 
beyond the reach of experiments in the foreseeable future [2, 3]. Thus if the 
forthcoming [3] or proposed experiments detect a nonzero d^, it will be an 
indication of physics beyond the SM. 

Recent neutrino data proves that Lepton Flavor (LF) has been violated in 
nature. The effect can be explained by mixing in the neutrino mass matrix. 
In principle such a LF Violation (LFV) in neutrino mass matrix can lead to 
the LFV decays, r — >■ 67, r — > /i7 and fi ^ e'-f [4]. However, if the neutrino 
mass matrix is the only source of LFV, the branching ratio of these decays 
will be so small that cannot be detected in foreseeable future. In the future, 
if experiments report a nonzero branching ratio for any of the aforementioned 
LFV decays [5] , we will conclude that the SM has to be augmented to include 
more sources of LFV. In the context of Minimal Supersymmetric Standard 
Model (MSSM), which is arguably the most popular extension of the SM, 
there are several sources for CP- violation as well as for LFV which can lead to 
effects exceeding the present experimental bounds. The bounds on Bic{ij 
ii'j) constrain the sources of LFV in the MSSM. Moreover, the bounds on 
the EDM of the elementary particles constrain the CP-violating phases of 
the MSSM. For vanishing LFV sources, the bounds from the EDMs on the 
CP- violating phases of MSSM parameters have been extensively studied in 
the literature (for an incomplete list see) [6, 7, 8]. In [9], the effects of 
the phases of LFV masses as well as the LFV trilinear ^-couplings on 
have been studied. However, [9] does not discuss the possible effects of the 
phase of A-,- (the trilinear supersymmetry breaking coupling of stau) . Notice 
that although A^- is a LF conserving couphng which deals only with the 
staus, in the presence of LFV, it can affect the properties of leptons of other 
generations. 

In this paper, taking into account the possibility of the LFV in soft super- 
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symmetry breaking terms, we will focus on the possible effects of the phase 
of Aj. on the electric dipole moment of the electron. As is well-known, the 
phase of Ar, (pA-rt can also manifest itself in the decay and production of 
staus [10]. One of the goals of the proposed state-of-the-art ILC project is 
detecting such effects [11]. It is therefore very exciting to learn about the 
value of (pA-r by present or forthcoming low energy experiments. 

We show that for Br(r 67) close to its present bound, the bound on 
de can already constrain . We discuss how other sets of the CP- violating 
phases can mimic the effect of on dg and suggest some ways to resolve 
the degeneracies. Recently, it has been shown in [12] that by measuring the 
spin of the final particles in the LFV rare decays, one can extract information 
on the CP-violating phases of the underlying theory. In the present paper, 
we however do not take into account such a possibility and focus on the 
spin-averaged decays rate. 

The paper is organized as follows. In Sec. 2, we specify the model that we 
study in this paper and summarize the observable effects that can be used to 
extract information on the parameters of the model. In Sec. 3, we discuss how 
CP- violating and LFV parameters affect de and other observable quantities 
and present scatter plots that explore the parameter space. In Sec. 4, we 
first enumerate the possible CP- violating phases and evaluate their respective 
effects with special emphasis on the possibility of cancelation. Section 5, 
summarizes our conclusions. The formulae for calculating the rate of LFV 
rare decays and d^ are summarized in the appendix. 

2 The model and its observable effects 

In this section we specify the model and the sources of LFV and CP- violation 
that we are going to study. 

In this paper, we consider the Minimal Supersymmetric Standard Model 
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with superpotential 

W^MSSM = -Yie^i Li - Hd- II Hu- Hd (1) 

where Li, Hu and Hd are doublets of chiral superfields respectively associated 
with left-handed leptons and the two Higgs doublets of the MSSM. In the 
above formula, e^^ is the chiral superfields associated with the right-handed 
charged lepton fields. The index "i" determines the flavor; i = 1, 2, 3 respec- 
tively correspond to e, /i, r. Notice that we have chosen the mass basis for 
the charged leptons [i.e., Yukawa coupling of the charged leptons is taken to 
be diagonal). Notice that the Yukawa terms involving the quark supermul- 
tiplets have to be added to (1). However, since we are going to concentrate 
on the lepton sector, we have omitted such terms. At the electroweak scale, 
the soft supersymmetry breaking part of the Lagrangian in general has the 
form 

L^^^ssM ^ _ (^MiBB + M2WW + H.c.) 

- {AiYi5ij + Aij)e%- Lj ■ Hd + H.c.) - (mlJijLj - e^J {mljije'iij 

- mj^^HlH^- ml^HlHd-{BH H^-Hd + R.c.), (2) 

where the "i" and "j" indices determine the flavor and Lj consists of (i>j en)- 
Notice that we have divided the trilinear coupling to a flavor diagonal part 
(AiYiSij) and a LFV part {Aij with An = 0). Again terms involving the 
squarks as well as the gluino mass term have to be added to (2). The Her- 
miticity of the Lagrangian implies that m^^ , and the diagonal elements 
of and are all real. Moreover, without loss of generahty we can 
rephase the fields to make the parameters M2, Bh as well as Yi real. In such 
a basis, the rest of above parameters can in general be complex and can be 
considered as sources of CP-violation giving contributions to EDMs. 

After electroweak symmetry breaking, the A-terms in (2) as well as the 
terms in superpotential induce left- right mixing. The Hermitian 6x6 mass 
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matrix of (e/j)j and (6^)^ can in general be written in terms of three 3x3 
mass sub-matrices m|, mj^ and as follows 

^Dd)- (3) 

The formulae for m|, m|, and in terms of the soft supersymmetry 

breaking potential are given in Eqs. (12,13,14) of the appendix. With above 
Lagrangian and superpotential, the LF is conserved if and only if Aij = 
and the off-diagonal elements of m?^ and m?^ vanish (for i ^ j, {m'j^)ij = 
{m'j^)ij = Aij = 0). At the one loop level, in the lepton conserving case, 
each A term can contribute to the EDM of only the corresponding fermion. 
For example, at the one loop level, the phase of Ar will have no effect on 
df. but can induce an EDM for the tau lepton of order of Im(A^)m^/mgugY- 
Considering the fact that the bound on the EDM of tau is much weaker than 
this [13], no bound on (f)A^ from dr can be derived. In the LF conserving 
case, de will receive significant contributions from the phases of Ae, fi and 
Ml. Thus, the strong bound on d^ can be translated into bounds on the 
phases of these parameters. In the literature, it is shown that for relatively 
low scale supersymmetry (msusY < 500 GeV), the bounds on these phases 
from de are very strong [8] unless severe cancelation takes place [6]. 

At the two-loop level, even in the lepton flavor conserving case, the phase 
of At- can induce a contribution to dg as well as to dn [14]^. For relatively 
large values of tan P (tan P > 10) and msusY — 100 GeV, the bound on dg 
can be translated into a bound of order of few hundred GeV on ImfA,-]. The 
limit from the bound on dn even is less stringent ^. 

For LFV case, the A-term associated with a definite lepton flavor can in 
principle affect the EDM of a lepton of another flavor, even at the one-loop 

^Although in [14] the two-loop effects of only Im[Ab] and Ini[At] on dg and rf„ have 
been discussed, similar discussion also holds for Im[A7-]. 

^We would like to thank the anonymous referee for pointing out such a possibility. 
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level. In particular if the er element of m| and m|. or Aij are nonzero, the 
phase of Aj. can induce an EDM for the electron exceeding the present bound 
by several orders of magnitudes. As a result, in this case the strong bound on 
dg can severely restrict the phase of At-. In order to study this bound, we have 
to first consider the bounds on the LFV masses and A-terms from the bounds 
on the LFV decay modes of the charged leptons. Notice that throughout this 
paper we have implicitly assumed that the origin of LFV lies at an energy 
scale far above the scale of the electroweak symmetry breaking. We therefore 
have the same LFV- violating elements for the left-handed charged lepton and 
sneutrino mass matrices. 

The strongest upper bound on the LFV elements of the slepton mass 
matrices comes from the following experimental bound: 

Br(/x ^ 67) < 1.2 X 10"^^ (4) 

which for msusY ~ 100 GeV implies (mDe;,, (m|)e;. ~ 10"^ - lO-^mlysY) 
and Ae^,A^e ~ 0-05?7i|(jgY/(i?d). Throughout this paper we will set these 
LFV elements equal to zero: 

{ml)eiu = {mR)eu. = and = = 0. 

There are also strong bounds on the branching ratios of LFV decay modes 
of the tau lepton: 

Br(r ^ 67) < 9.4 x 10"* [15] (5) 

and 

Br(T ^ //7) < 1.6 X 10-^ [15] (6) 

which can be respectively translated into bounds on the re- and TyU-elements. 
However, it can be shown that the bound on the re-elements from (5) are 
not very strong and these elements can be of the same order as the diagonal 
elements. Suppose that both the re- and r/x-elements are nonzero. This 
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means the e- and /x-lepton numbers are both violated and the /i 67 decay 
can therefore take place despite the vanishing ^e-elements. In fact, for rela- 
tively large re-elements saturating the bound from (5), the bound (4) can be 
translated into a strong bound on the r/^-elements which is more stringent 
than the bound from Br(T /ij). Throughout this paper, we will set all the 
T/i equal to zero: 

{iT^DnT = {'m%)ij,r = and = = 0. 

In this scenario, the yU-fiavor number is thus conserved. 

As shown in the literature, integrating out the heavy supersymmetric 
states, T —>■ ej can be described by the following effective Lagrangian 

eeimrq/3 [eRa''^{AL)erTL + eL(7"^(^i^)erTfl] + H.c. (7) 

where is the photon field and qp is its four-momentum and cr"^ = |[7", 7^]. 
The formulae for Al and Ar in terms of the supersymmetric parameters 
have been derived in [16] for the CP-invariant case. We have re-derived the 
formulae for the CP- violating case. The results can be found in the appendix. 
Our results are in agreement with [16] in the CP-invariant limit. Using (7) 
it is straightforward to show that in the rest frame of the tau lepton, the 
partial decay rate is given by 

'^^^IZr^ = [(l(^^)e.r(l + cos^) + \{An)Mi - cos9)] (8) 

where 9 is the angle between the spin of the tau and the momentum of the 
emitted electron. Integrating over 6 we obtain 

r(r - e7) = ^mU\{A,W + \{AnW). (9) 

Notice that different sets of LFV mass matrix elements can result in the same 
rate for r — > 67. Let us define the Ap asymmetry as follows 

fn '^^i' ^o'^ dcosO — r°-, d COS 

4 X dcosU J —1 acosO flO") 

r{r — >• 67) 



7 



Using (8), we can prove that 
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Thus by measuring the partial decay rate of r, we will be able to extract not 
only (|(>lL)er|' + \{AR)er?) but also (l(AL)erl' " \{AR)er?)- In [17], it has 
been shown that by studying the angular distributions of the final particles 
at an e~e+ coUider such as a B-factory, it will be possible to derive Ap. Ap 
provides us with more information on the LFV parameters of the underlying 
theory. For example, if the source of LFV is a canonical seesaw mechanism 
embedded in the MSSM, we expect (m|j)er ^ {nT'DeT and A^r <^ A-j-e which 
means |(^L)eTp <^ |(^ii)erP and therefore Ap — > —1. In this paper, we 
will study the correlation between de, Br(T — > cy) and Ap and discuss the 
possibility of resolving the degeneracies by combining the information on 
their values. 



3 New contributions to de in the presence of 
LFV 

Let us for the moment suppose A^.^ — A^^ — 0. As illustrated in Fig. 1, for 
nonzero {'ni]^)eT and {m\)re-i the phase of At can induce a contribution to de- 
As a result for definite values of the off-diagonal mass elements, the bound on 
de can be interpreted as a bound on (j)A^ or on Im^Ar). Consider the case that 
both {mp)Te and {m\)re are close to the corresponding bounds from Br(r — > 
67). In this case, Br(r — > 67) is close to its present bound and Ap takes a 
value in the interval (-1,1); i.e., Ap ^ ±1. For such a configuration, we expect 
the bound on Im(AT-) to be more stringent than the bound on Im(Ae) because 
the effect of Im(A^) is given by m^Im(A^)(m^)eT(^i)er/'"^Lsj;, whereas the 
effect of Im(Ae) is proportional to mJxa{Ae) / rril^gy. Now, suppose only one 
of {rn^pjeT and {m\)Te is close to its upper bound from Br(r 67) and the 
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other is zero or very small. In this case, Ap will converge either to 1 (for 
(m|j)eT- ^ (m|)eT-) or to -1 (for {m\)er ^ ("^R)er)- From Fig. 1 we observe 
that if only one of {m^^er or {m'j^)re is nonzero and the rest of the er entrees 
(including A^r and A^g) vanish, at one-loop level, the phase of A^- cannot 
contribute to dp. 




Figure 1: A neutralino exchange diagram contributing to de- The photon 
can attach to any of the e^, f^, or propagators. The boxes on the 
left and right sides respectively depict insertion of {m\)er and (m|.)eT- The 
circles indicate insertion of the A^- vertex and the vacuum expectation value 
of if,. 



Figs. 2-4 demonstrate this observation. To draw the figures in this pa- 
per, we have chosen the mass spectrum corresponding to the a benchmark 
proposed in [18]. However, we have allowed the mass spectrum of the staus 
to slightly deviate from these benchmarks. Notice that at these benchmarks, 
the lightest stau is considerably heavier than the lightest neutralino so stau- 
neutralino coannihilation cannot play any significant role in fixing the dark 
matter relic density. As a result, a slight change of stau parameters will not 
dramatically affect the cosmological predictions. Although for illustrative 
purposes we have displayed the mass insertion approximation in Fig. 1, to 
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calculate de and Br(r 67) we have used the exact formulae (without the 
mass insertion approximation) presented in the appendix. 

Fig. 2 depicts rfg in terms of the sine of (pA^ for Aij = and various 
values of {m\)eT and (m|,)eT-- This diagram demonstrates that for {m\)eT and 
(m|j)eT close to their bounds from Br(T 67), a very strong bound on (p^-r 
can be derived. That is while if there is a hierarchy between these elements, 
the bound will be much weaker. Notice that for the input parameters chosen 
for this figure, Br(T 67) lies close to its present bound: 10~^ < Br(T — > 
67) < 10-^. 

Fig. 3 demonstrates the correlation between Ap and d^.. As explained in 
the caption, the input mass spectrum is that of the a benchmark [18] and 
\Aj-\ — 500 GeV. We have set A^t and A^-e equal to zero and the maximal 
value for the CP-violating phase is chosen: = 7r/2. (m|),-e and (m|j)^e 
pick up random values at a logarithmic scales. Points for which Br(T — > 
67) exceeds its present bound are eliminated. Fig. (3-a) shows us that if 
(pAr — 7r/2, de for a significant portion of the scatter points exceeds the 
present bound. Fig. (3-b) shows Ap versus Br(T — > 67) for the same scatter 
points. To illustrate the correlation between Ap and dg, we have shown the 
corresponding scatter points in Fig 3-a and Fig 3-b with the same color and 
symbol. That is at points marked with green dots, de exceeds its present 
bound and at the scatter points marked with blue cross "x" 10~^^ < de < 
1.4 X 10~^^ e cm. The scatter points depicted by pink circle, which appear 
in Fig. (3-b) as two pink horizontal lines at Ap — ±1, correspond to de < 
10~^^ e cm. From Fig. (3-b) we conclude that for Aer — A^-e — 0, the bound 
on de can be satisfied if either Br(T — > 67) is very small (which means that all 
the LFV masses are very small) or is close to ±1 (meaning that there is 
a hierarchy between the LFV elements). In other words within this scenario, 
if future searches find 5 x 10~^° < Br(T 67) and -0.9 < Ap < 0.9, 
the bound on dg should be interpreted either as a bound on (pA^ or as an 
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Figure 2: de versus sin (/)^^. The input parameters correspond to the a bench- 
mark proposed in [18]: = 375 GeV, mo = 210 GeV, M1/2 = 285 GeV 
and tan/? = 10 and we have set |At-|=500 GeV. All the LFV elements 
of the slepton mass matrix are set to zero except {mj^)eT and (m^)er- 
The dotted (pink) line labeled (a) corresponds to (m|)er=3500 GeV^ and 
(m|j)er=15000 GeV^. The dashed (green) line labeled (b) corresponds to 
(mi)er=50 GeV^ and (m|)er=37000 GeVl The sohd (red) line labeled (c) 
corresponds to (m|)e^=3500 GeV^ and (m^)er=30 GeV^. The horizontal 
doted line at 1.4 x 10~^^ e cm depicts the present experimental limit [13] on 

de- 
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Figure 3: a) Scatter plot of de versus Br(r 67). The input parameters 
correspond to the a benchmark proposed in [18]: = 375 GeV, mo = 
210 GeV, Ml/2 = 285 GeV and tan/3 = 10. We have however set = 7r/2 
and l^^l = 500 GeV. All the LFV elements of the slepton mass matrix are 
set to zero except (mDe,- and (m|,)eT which pick up random values at a 
logarithmic scale respectively from (5.9 x 10~^ GeV^,5.9 x 10^ GeV^) and 
(3.7 X 10-3 GeV^3.7 x 10^ GeV^). The horizontal line at 1.4 x lO'^^ e cm 
depicts the present experimental limit [13] and the one at 10^^^ e cm shows 
the limit that can be probed in the near future [2]. b) Scatter plot of 
versus Br(r — > 67). For each scatter point in Fig. 3-a there is a counterpart 
in Fig. 3-b corresponding to the same input values for the er elements which 
is shown with the same color and symbol. Notice that points shown in pink 
(corresponding to < 10"^^ e cm) all lie on the Ap=±l. 
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Figure 4: Similar to Fig. 3 except that (m|)er = (^fe)eT = and instead 
(m|^)er(= AeriHa)) and {m'j^j^)re{= A^e{Hd)) pick up random values at a 
logarithmic scale from (1.2 x 10"^ GeV^, 1.2 x 10^ GeV^). For each scatter 
point in Fig. 4-a there is a counterpart in Fig. 4-b corresponding to the same 
input values for the er elements which is shown with the same color and 
symbol. 

indication for a cancelation between different contributions from (pA-r and 
other possible CP-violating phases. In the next section, we shall elaborate 
on the latter possibility in more detail. We emphasize that to draw this 
conclusion we have taken Aer = Are = 0. Repeating the same analysis with 
different benchmarks we have found that the pattern and the correlation 
among Br(r — 67), Ap and de are not sensitive to the input values of LF 
conserving parameters. 

To draw Fig. 4, {m\)eT and {m'j^)er are set equal to zero and instead 
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Figure 5: Similar to Fig. 3 except that here in addition to {m\)f.^ and (m^)er, 
(m|^)eT-(= ^er(-f^(i)) and {jn\^)re{= ^reiHd)) are allowed to be nonzero. 
The values of {m\)eT and {fnj^)eT are randomly chosen respectively from 
(0.59 GeV^ 5.9 x 10^ GeV^) and (3.7 GeV^ 3.7 x 10^ GeV^) at a logarithmic 
scale. {jn\^eT and {jn\j^re pick up random values at a logarithmic scale 
from the interval (0.12 GeV^, 1.2 x 10^ GeV^). For each scatter point in 
Fig. 5-a there is a counterpart in Fig. 5-b corresponding to the same input 
values for the er elements which is shown with the same color and symbol. 
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random values for Aer and Are are taken. The LF conserving parameters 
for Figs. (3) and (4) are the same. Like the case of Fig. (3), a significant 
portion of the scatter points have exceeding the present bound. Moreover 
a similar correlation between dg, Ap and Br(r — > 67) emerges. That is points 
marked with green dot (corresponding to rfg >1.4 x 10"^'' e cm), with blue 
"x" (corresponding to 10^^^ < dg. < 1.4 x 10^^'^ e cm) and pink circles 
(corresponding to de < 10^^^ e cm) are scattered respectively from right to 
left. Notice however that in contrast to Fig. (3-b), Fig. (4-b) includes scatter 
points with —0.9 < Ap < 0.9 and Br(r — > 67) ~ 10~^ that satisfy the present 
bound on d^ (the points marked with "x" in the plot). We have repeated 
the same analysis with the S benchmark and have found a similar pattern. 
Similarity between the patterns means that the above observation does not 
depend on the input values for the LF conserving parameters. 

In Fig. 5, the input values for the LF conserving parameters are taken 
to be the same as those for Figs. 3 and 4. However, (m|)er, ("T.|j)er, ^er 
and Are all take nonzero random values. Fig. (5-a) contains features of 
both Figs. (3) and (4). The significant point is that setting all the er mass 
elements nonzero, the correlation among Ap, dg and Br(T 67) is lost. That 
is Fig. (5-b) contains points with Br(T — > 67) ~ 10~^, —0.9 < Ap < 0.9 and 
de < 10~^^ e cm. The presence of these points can be explained by the fact 
that when Aer and {m\)eT are nonzero but A^-e — {frtpjeT — (or equivalently, 
when A^e and {m?p)eT 7^ but Aer — {'mj^)eT — 0) de, despite large 0^^, 
remains zero but Al can be of order of Ar which yields —0.9 < Ap < 0.9. 
As a result, without independent knowledge of the ratios of LFV elements, we 
cannot derive any conclusive bound on even if we find —0.9 < Ap < 0.9 
and 10~^ < Br(T — > 67). We have repeated the same analysis for other 
benchmarks and the results seem to be robust against changing the mass 
spectrum. 

As explained earlier, some models predict a certain pattern for LFV. For 
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example, within the framework of the seesaw mechanism embedded in the 
constrained MSSM, we expect the LFV to be induced mainly on the left- 
handed sector [19]. That is we expect (w^/j)er ^ ("^^i)eT and Agr/^re ^ 
rrie/mT- <C 1. This model predicts Ap = —1. On the contrary, within the 
supersymmetric SU(5) GUT model without right-handed neutrinos, the LFV 
is induced only on the right-handed sector [20] which implies Ap — 1. For 
both of these cases, de induced by is negligible. 

In the above discussion, we have used the bound on (f)A^ and on Im(A,-) 
interchangeably. To clarify the relation between these two. Fig. (6) has been 
presented which shows curves of — lA x 10"^^ e cm (the present bound) 
for the a and S benchmarks and various values of the off-diagonal elements. 
The values of the LFV elements are chosen in a range to obtain Br(r — > 67) 
close to the present bound; i.e., 10"^ < Br(r — > 67) < 10~^. Each curve 
can be considered as the upper bound on 0^^. These figures also confirm 
that when there is a hierarchy between the left and right LFV elements, the 
bounds are weaker. As expected, the curves have a shape close ImlA-,-) — 
\At \ sin (f)A.r — cte. 

In summary, within a model that A^j — 0, if Br(r 67) turns out to be 
close to its present bound and Ap deviates from +1 and -1, the bound on dg 
puts a strong bound on Im(A,-). However, if Ap — ±1, the bound on dg can 
be explained by a hierarchy between the (mD^-e and {mj^jre elements instead 
of by the smallness of Im(^^). Similar discussion holds for the scenario in 
which (m^)er = ("^i)eT — and instead A^t- and A^^ are nonzero: while for 
Ap 7^ ±1, the bound on de can severely restrict (pA^, for = ±1 we cannot 
obtain any bound on from dg. However, within a scenario that Aer, Are, 
{mj^jre and {rri^)re are all large, we cannot derive any bound on even if 
Ap ^ ±1. Thus, in order to derive a conclusive bound on 0^^, one has to 
resolve these degeneracies seeking help from an experiment other than the 
rare r decay. Studying LFV signals at a e~e'^ collider with energy of center 
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Figure 6: Contour plots for de=lA x 10 a) The a benchmark pro- 
posed in [18] is taken as the input. The dotted blue curve corresponds to 
(mi)er=3500 GeV^ (m|)er=1500 GeV^ and (miR)^e=(mi^)e^=0 and the 
thick solid red curve corresponds to (m|)er=3000 GeV^, (m|.)er=1000 GeV^, 
(m|^)er=300 GeV^ and (m|^)T-e=100 GeV^. The dashed green curve cor- 
responds to (m|)e^=3000 GeV^ (m|)e^=100 GeV^ (m|^)e^=10 GeV^ and 
(m|^)^e=500 GeV^. b) The 6 benchmark proposed in [18] is taken as the 
input: = 920 GeV, mo = 500 GeV, M1/2 = 750 GeV and tan/3 = 10. 
The dotted blue curve corresponds to (m|)e^=7 x lO'' GeV^ and (m|j)e^=2 x 
10^ GeV^ and the thick red curve corresponds to (m|)er=2 x 10^ GeV^, 
(m|),^=3 X 10^ GeV^ {ml^)er=SOOO GeV^ and (m2^)^e=7000 GeV^ The 
dashed green curve corresponds to (m|)er=lx 10^ GeV^, (m|j,)er=3000 GeV^, 
(^Der=30 GeV^ and (mi^)^e=8000 GeV^ 
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of mass of a few hundred GeV can help in this direction [21]. In this paper 
we have concentrated on the possibihties that the ongoing experiments can 
bring about. Studying the possibihties with ILC is beyond the scope of the 
present paper. 

4 Degeneracies between different sources of 
CP-violation 

In the previous section, we had assumed that the only source of CP-violation 
is the phase of A,,-. However, within the framework of general MSSM, there 
are multiple sources of CP-violation. In the basis described in Sec. 2, these 
phases include the phases of Ae, the //-term and Mi (the Bino mass) that can 
contribute to de regardless of the conservation or violation of LF. Within the 
scenario considered in this paper, in addition to these sources, the phases of 
{ml)eT and (m|j)er can be also considered as independent sources 
of CP- violation that can contribute to de- If we assume that only one of 
these various phases is nonzero, the present bound on de can be interpreted 
as a strong bound on the nonzero phase. However, in general when more 
than one phase is present, the effects of different phases can cancel each 
other [6]. Moreover, if the forthcoming searches report a nonzero de, without 
additional information, we cannot disentangle the source of CP- violation. In 
this section, we discuss the degeneracies between different sources of CP- 
violation with special emphasis on the possibility of cancelation. 

Figs. (7-9) display the degeneracies between possible CP- violating phases. 
To draw these figures, we have inserted the mass spectrum of the a bench- 
mark proposed in [18] (see caption of Fig. 3 for the values of the relevant 
parameters) and we have set A^ — Ae — 500 GeV. Each of Figs. (7-9) cor- 
responds to a different set of absolute values for the LFV elements (m|)eT, 
(m|j)eT, Are and Aer- Each curve in these figures shows de versus Br(T — > 67) 



18 




0.2 



0.4 



0.6 



0.8 



Br(T^eY) [10" 



Figure 7: df. versus Br(r 67) as the CP-violating phases vary between zero 
and TT. The input parameters correspond to the a benchmark proposed in 
[18]:|/i| = 375 GeV, uiq = 210 GeV, M1/2 = 285 GeV and tan/? = 10 and we 
have set y4^=y4e=500 GeV. All the LFV elements of the slepton mass matrix 
are set zero except that | (m|)er |=3500 GeV^ and |(m|,)er| = 15000 GeV^. To 
draw the curves all phases are set zero except one that varies between and n. 
As illustrated in the legend of the figure, the thin solid red curve, dotted grey 
curve and light blue dashed curve respectively correspond to varying phases 
oi At-, h and Mi. The thin black and thick pink dotted curves correspond 
to the phases of {m\)eT and (m|,)e^ which for A^.^. = A^-e lie over each other. 
The thick green vertical line stretching up to de=2 x 10~^^ e cm depicts the 
effect of the phase of A^. The horizontal line at 1.4 x 10"^'' e cm depicts 
the present experimental limit [13] and the vertical line shows the present 
experimental bound on Br(r 67) at 9.4 x 10~^ [15]. 
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as a certain CP-violating phase varies from zero to vr while the rest of the 
phases are set to zero. As expected all the curves converge at rfg = which 
corresponds to the zero value of the varying phase. As the value of the vary- 
ing phase reaches 7r/2, obtains its maximum value so the peak of each 
curve corresponds to the varying phase equal to 7r/2. The horizontal lines 
at de = 1.4 X 10~^^ e cm in the figures show the present upper bound on d^ 
and the vertical lines at Br(r — > 67) = 9.4 x 10~^ show the present bound 
on Br(r 67). In the following, we discuss these figures one by one. 

Drawing Fig. 7, we have set |(m|)e^| = 3500 GeV^ |(m|)er| = 15000 GeV^ 
and Ag.T = ^re = 0. The CP-violating phases that can contribute to dg in- 
clude (f)Ae, <t>Ar-i (pfj., 0Mi and the phases of {m\)eT and {m^^eT- The thick 
vertical fine at Br(T 67) = 3.8 x 10~^ corresponds to the variation of 
in [0,7r]. This line shows that Br(r — > 67) does not significantly change 
as (t)Ae varies. The reason is that the effect of on Br(T — > 67) is much 
smaller than the dominant effect. The line associated with (the thick 
line) reaches values of de up to one order of magnitude higher than the present 
bound on d^ which means if is the only contributor to de, it cannot be 
larger than O(O.l). This bound is similar to the bound in the LF conserving 
case. Notice that the effects of the rest of phases can exceed the maximal 
contribution from by more than one order of magnitude. In this figure, 
the curves associated with the phases of (m|)er and (m|j)er, which are de- 
picted by black and pink dotted curves, coincide. This observation is valid as 
long as — ^re — because the diagram shown in Fig. 1 - which in this 
case is the only diagram contributing to d^ - is sensitive only to the relative 
phase of {m\)eT and {m^^er- Another pecufiar feature of Fig. (7) is that 
the contribution of to de can exceed the maximum de from 0^. This is 
opposite to the LF conserving case in which the effect of 0^ is larger because, 
while 0Mi can induce de only through the subdominant neutralino-exchange 
diagram, 0^ can induce EDM also through the dominant chargino-exchange 
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diagram. In contrast to the LF conserving case, in the case of Fig. (7) the 
neutrahno exchange diagram dominates because as explained in the previous 
section, once we turn on the er elements, the neutralino-exchange diagram 
contributing to rfg is enhanced by a factor of mT-/me. As a result, the effect 
of (pMi is enhanced. 

Now let us discuss the degeneracy and the possibility of cancelation among 
different contributions. Replacing a phase with its opposite value, its con- 
tribution to (ie will change sign. As a result if we find two phases whose 
contributions to \de\ have the same values, we can conclude that cancelation 
can take place for at least one pair of values. Fig. 7 shows that the curve 
associated with (pAe has a complete overlap with the low phase part of the 
other curves. That is for any value of (j)A^, there is a value for other phases 
which can mimic the effect of (j)A^- Thus, if the future EDM searches re- 
port a nonzero value for de, there will be an ambiguity in interpretation of 
the observation in terms of the phases. Cancelation is another consequence 
of this overlap. This figure shows that turning on more than one nonzero 
phase, cancelation can make even the maximal value of consistent with 
the present bound on d^. That is the contribution of — 7r/2 can be can- 
celed out by the effect(s) of any of the phases 0Mi, 0;ti or the phases of (m|)er 
and (m|j)er if these phases are ~ 7r/500. The contribution of = ±7r/2 
can be also canceled out by the contribution of (f)Ar if \4>Ar \ — 7r/80. 

Whereas the phase of Ag can only show up in d^, the rest of phases can 
manifest themselves as CP-odd effects at ILC [11]. Moreover, and 0^ 
can give a detectable contribution to the EDM of the neutron [22] , mercury 
[23] and deuteron [24] through inducing chromoelectric dipole moments and 
EDMs to the light quarks. Thus, from the experimental point of view, can- 
celation between the effects of these phases is more exciting as it can open 
up the possibility of large phases and therefore CP-odd observable quantities 
in experiments other than de searches. Fig. 7 shows that there are values 
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of 0Mi and/or the phases of (m|)er and/or {iTij^)er whose contribution to 
df. can cancel even the maximal effect from 0^. That is even 0^ = ±7r/2 is 
consistent with the bound on but in order to cancel the effects down to the 
present bound a fine-tuning better than 0.1% is needed. In [9], it was also 
shown that turning on the LFV elements of mass matrix, the effect of 0^ on 
de can be canceled by the effects of the phases of the off-diagonal elements. 
This result is obviously sensitive to the largeness of the absolute values of the 
LFV elements which are constrained by the null results of searches for the 
LFV rare lepton decays. Bearing in mind that since [9], these bounds have 
significantly improved the above discussion can be considered as an update 
and re-confirmation of the claim in [9] in view of the recent bounds. Notice 
however that a large 0^ can also give contributions to d„ and dng exceeding 
their present bounds. The correlation between dg and dug has been system- 
atically studied in [25]. To satisfy the bounds on dn and dng in the presence 
of a large (p^, there should be another cancelation scenario at work in the 
quark sector. This further suppresses the allowed parameter space; i.e., a 
double-folded fine-tuning. 

Let us now suppose that there is a symmetry or a mechanism that sets 
— 0Mi = so the bounds on dn and dng are naturally satisfied without the 
above mentioned double-folded fine-tuning problem. Let us also suppose that 
(f)Ar and the phases of (m|)er and (m|j)eT are large. Fig. 7 shows that this 
scenario is not ruled out by the d^ bound because there is still the possibility 
of cancelation between the contributions of the nonzero phases. To cancel 
the effects of (pA^ — 7r/2 on de down to the present bound, a fine-tuning 
better than 1% is required. 

Fig. (8) has an input similar to that of Fig. (7) except that A^t and A-j-e 
are set nonzero and smaller values for (m|)eT and (m|j)eT are chosen. Notice 
that unlike Fig. (7) curves associated with the phases of (m|)er and (m^)eT 
split. The peak of the (pA^ curve in Fig. (8) lies one order of magnitude below 
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Figure 8: Similarly to Fig. 7 except that here we have set |("^i)eT| = 
1000 GeV^ |(m|)e^|=5000 GeV^ and \{m\^^r\ = K^Drel = 300 GeVl 
The thin solid red curve, light dash-dotted grey curve, light blue dashed 
curve, sohd dashed dark blue curve and thick black dotted curve respectively 
correspond to the varying phase of Ar, /i. Mi, (m|)er and (m|.)er- The 
thick green vertical line stretching up to de=2.1 x 10~^^ e cm depicts the 
effect of the phase of Af,. The light pink curve and thin green curve respec- 
tively correspond to phases of (^|/j)Te and {m^^eT- The horizontal line at 
1.4 X 10^^^ e cm depicts the present experimental limit [13] and the vertical 
line at 9.4 x 10~^ shows the present experimental bound on Br(r 67) [15]. 
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that in Fig. (7). That is because the absolute values of (m|)er and (m|.)eT 
in Fig. (8) are smaller. Had we set these elements larger, the effect of 
on de would have been larger but also the value of Br(r ej) would have 
increased. The rest of the argument for Fig. (7) holds for Fig. (8), too. 

Fig. (9) displays the dependence of and Br(r — > 67) on different phases 
for the case that there is a hierarchy between the left and right LFV elements: 
\Aer\ -C lArel and \{m'j^))eT\ <^ |('^i)er|- Bccausc of this hierarchy, the effect 
oi(l)Ar on de has dropped below the present bound which is expected following 
the discussion in the previous section. The lines associated with the phases 
of {m'j^)eT and {m\p()eT appear as vertical lines which means Br(T — >• 67) 
does not depend on these phases. This is expected because |(m|j)er| and 
|(m|^)e^| are very small. However, the effects of their phases can still exceed 
the present bounds. The figure also shows that Br(T — > 67) strongly depends 
on the phases of (m|)er and {m\j^re- The effect of is similar to the 
previous cases. The 0Mi curve also appears as a vertical line which means 
Br(T — > 67) does not strongly depend on The effect of 0Mi on de in case 
of Fig. (9) is one order of magnitude smaller than the case of Fig. (8) and, like 
the LF conserving case, is smaller than the effect of 0^. In contrast to Figs. 
(7) and (8), in this case 0^ ~ 7r/2 is ruled out by the bound on de because 
the effects of other phases will not be large enough to cancel the effect of 
0^ ~ 7r/2. However, for 0^ < 7r/30 the effect of 0Mi and for 0^ < 7r/500 the 
effects of the phases of LFV mass elements as well as that 0^^ can cancel the 
contribution from 0^ to de- 

5 Concluding remarks 

In this paper, we have discussed the effects of the phase of trilinear ^-coupling 
of the staus, 0^,., on de in the presence of nonzero LFV er elements of the 
slepton mass matrix. We have shown that for a large portion of the parameter 



24 



1000 - 



100 



E 
o 

CD 



T3 



10 



r 



1 - 



0.1 



0.01 



0.2 



0.4 



0.6 



0.8 



— I r 

^m^R mil 

(m2LR)xe 
(m LR)eT 





1.2 



Br(T^eY) [10"'] 



Figure 9: Similarly to Fig. 7 except that here |("^|)er| = 3000 GeV^, 
\im%)er\ = 50 GeV^ |(mi^)er|= \A,r\{Hd)=3 GeV^ and |(mi^).e| = 
\Are\{Hd)=4:00 GeV^. The thin solid red curve, light dotted grey curve, 
thin dotted blue curve and thin dotted pink curve respectively correspond 
to varying phases of Ar, fi, {m\)(,T and {jn\^re- The pink, green, black and 
dark blue thick vertical lines at Br(r — 67) = 7.5 x 10~^ (which reach up 
4 = 6.8 X 10-27,2 X 10-26,3.2 x 10-^^, 1.1 x lO^^s e cm) depict 4 versus 
Br(r — >• 67) as the phases of respectively {jn\^eT) A^, {rn'j^)er and Mi vary 
between and vr. The horizontal line at de = 1.4 x 10-^'' e cm and the vertical 
line at Br(r —>■ 67) = 9.4 x 10-^ show the present experimental limits [13, 15]. 
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space consistent with the present bound on Br(r 67), the contribution of 
(f)j^^ to de can exceed the present bound by several orders of magnitude. 
The effect of (f)A^ on strongly depends on the ratios of the LFV slepton 
masses (^i)er/(^R)eT and {m^Rler/i'i^'LRlTe- In other words, for a given 
Br(T 67) and (pAr = can take any value between zero and a 

maximum which depends on the value of Br(r — > 67) [see Figs. (3-a)-(5-a)]. 
We have shown that for specific case that {m\j^eT — {j^\R)re ~ [see Fig. (3- 
b)] or {m\)(,T = ij^RjcT — [see Fig. (4-b)], by measuring the asymmetry Ap 
defined in Eq. (10) we can solve this ambiguity. However, in the general case 
that all the er elements are nonzero, as shown in Fig. (5), the correlation 
between Ap and de is lost and to solve the ambiguity, extra information is 
needed. 

Assuming that is the only source of CP- violation contributing to d^ 
we have derived bounds on for various values of the LFV elements giving 
rise to Br(T — > 67) close to the present bound (see Figs. 2 and 6). We have 
then relaxed this assumption and discussed the possibility of cancelation 
between contributions of the different phases. We have shown that for large 
er mass elements saturating the present bounds, the effect of the phase of 
the Bino, on de is significantly enhanced which can be explained by the 
enhancement of the effect of the neutralino exchange diagram by a factor 
of mr/me- Taking into account the new bounds on branching ratios of the 
rare LFV tau decay, we have confirmed the results of [7] that with nonzero 
LFV effects cancelation scenario makes large values of 0^ consistent with 
the bound on de- We have discussed that the requirement to simultaneously 
satisfy the bounds on dg, dn and dug by cancelation imposes a double- folded 
fine tuning problem. 

We have shown that contributions from phases of (m|)er, {f^RjeTi ("^iR)eT 
and (w^|^)^-e can cancel the effect of on de- In summary, although in case 
of large er elements saturating the bounds from Br(r 67) , (pA^ can induce a 
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large contribution to rfg, still the possibility of cancelation and/or presence of 
a hierarchy between the LFV er mass matrix elements make even a maximal 
(pA-r consistent with the rfg bound even if Br(r 67) is found to be close to 
its present bound. Thus, still there is a hope to observe CP-odd effects at 
ILC [11]. 
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Appendix 

In this appendix, we summarize the formulas necessary for calculating Ap, 
Br(r — > 67) and d^. In this paper, we are interested in large LFV er elements. 
In this parameter range, the mass insertion approximation is not valid and 
one should work in the mass basis. Here, we first derive the coupling of the 
sleptons to neutralinos and charginos in the mass basis taking to account 
the CP-violating phases and mixing. We then present the formulas for Al 
and Ar defined in Eq. (7) as well as for the formula for de- Throughout this 
appendix we omit the spinorial indices for simplicity. 

In the flavor basis, the mass terms of (the superpartners of the left- 
handed charged leptons) and cr (the superpartners of the right-handed charged 
leptons) can be written as 

wo.=-(^i A)Mi{i)^-u 4)(4 ':|)(|) 

(11) 

where m| and m\ are 3x3 Hermitian matrices and is a general complex 
3x3 matrix. The elements of these matrices are as follows: 
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{ml)ij = {mljij + {ml)i6ij + m| cos 2/3(-- + sin^ Ow)Sij (12) 



= + ("^e)j<^u - "^1 cos 2/3 sin^ 6lvK(^y (13) 

and 

i'^LR)ij = nT'iiA - /U* tanP)6ij + Aij{Hd) (14) 

where m-^ and arc respectively the right-handed and left-handed slepton 
soft supersymmetry breaking mass matrices at the electroweak energy scale 
and Aij is the trilinear A-coupling [see Eq. (2)]. We can diagonalize the mass 
matrix of slepton by a 6 x 6 unitary matrix [/' as 

[U^Ml{UY%y = mld,y (15) 
The slepton mass eigenstate in terms of the chiral weak eigenstate are 

3 

e- = Efe^^^ + ^W«^] (16) 

i=l 

Since in the MSSM no ur exists, the neutrino mass matrix will be a 3 x 3 
matrix whose elements can be written as 

(mDij = {mljij + (^m| cos 2/3)5^,- (17) 

The mass eigenstate, Ux, is related to the weak eigenstate, uu, as 

3 

= 5^ t/.^ -^^^ (18) 

x=l 

Let us now consider the neutralino masses. The masses of neutralinos in the 
weak basis can be written as 

i^neutralino = -^(X°)^M^X° + H.C., (19) 
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where X^={B, W^, H^, H^) and 

(Ml —mzcpsw mzspsw \ 

M2 mzcpcw -mzSfsCw 

-mzCfsSw rnzc/scw -fx 

\ mzspsw -mzspcw -/i / 



(20) 



Here, sp — sin/3, c/3 — cos (3, sw — sin^^ and cw — cosOw- The mass 
matrix M^^ can be diagonahzed as follows: 

[ON*Mi^O^']AB = M^oJab (21) 

where Ojv is a unitary matrix and M^o are real positive mass eigenvalues. 
The mass eigenstates, x^, in terms of the weak eigenstates, X^, can be 
written as 

= {On)ab^. (22) 
The chargino mass terms can be written as 

i^chargino = -^(X±)^MgX± + H.C., (23) 



Mc-( n T ) (24) 



where {X^f = {W+ , H+, W', H') and 

c \ C Q 
with 

C^i J^' ^Sf^^wY (25) 

The chargino mass matrix C is a general complex matrix which can be diag- 
onahzed as 

W^CV"-' = dtag{\m~-\, |m~-|) (26) 
where and are unitary matrices that satisfy the following relations 

V%C^C)V'-' = dzag{\m~-\', \m^f) = U'{CC^)U'-\ (27) 
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Notice that we have defined If^ and in a way that the elements of the 
diagonal matrix W^CV^"^ are real positive. Eqs. (26,27) are invariant under 

^'^( ^'^( (28) 

Thus, there is an ambiguity in the definition of C/'^ and V'^ but the final results 
do not depend on the unphysical phases cti and a2, as expected. 
The mass eigenstates are related to the gauge eigenstates through 

Within the framework of the MSSM, the Icpton-slepton-neutralino cou- 
pling in the mass basis and in the Weyl representation can be written as 

6 

^2 = E ^uiN^Dxle. + e^NDx^Ae. + H.c, (30) 

x=l 

where the couplings are 

Y Z Iilw COb p 

NL = -^[2(Oiv):iitan^^t/^^^3 + -^^(0;v)l3^a- (31) 
The lepton-slepton-chargino coupling can be written as 

3 

= E ^iiiCL)xl% + eUctAjXA^.+a-C. (32) 

x=l 

where the couplings are 

CiAx = —fK ^ (33) 

y Imw cos p 

Let us now discuss the formulas for {^Ai)ij and (>lii)jj defined as 

eetm^eia"'5g^[(Ai)i,Pi + (AR)i,P^]e,- + H.c. (34) 

30 



where a"^ = |[7",7^] and qp is the four-momentum of the photon. Pl and 
Pr are respectively the left and right projection matrices. Notice that {Aijer 
and [A^jer defined in (7) are the er component of the 3x3 matrices {A^jij 
and {AR)ij. For the CP-conserving case the formulas for A^ and A^ have 
been developed in [16]. We have rederived the formulae for the CP- violating 
case. It is convenient to decompose Al and Ar as follows 

= Ai% + A^^^^ (35) 

where and A^^^^ respectively come from neutralino-slepton and chargino- 
sneutrino loops. In terms of the coupling in the mass basis we can write 

^^11 1 

X (1 - 6yA. + 3y\, + 2yL - %L In Va.) (36) 

MyO 1 

+ Nt,.NX-^T, ^(1 - Z/L + 2t/A.lnt/A.)], 

(1 - VAx) 

where uax = Af|o / and 

^ 3 1 1 I 



^^^^^^ 327r^ml' ^^^6{1-zax)' 

X {2 + Szax - + 'ozax In zax) 



m^i (1 - ZAxY 



+ C£i.C'f;,Tr^^:^ ^(-3 + ^ZAx - 4. - 21nz^.)] 



(37) 



ZAx^Ml_Jm\. Finally 



Notice that the forms of the above formulas are similar to those in [16]; 
however, the couplings A^^'-^ and C^'^ are slightly different because of the 
nonzero CP-violating phases. 
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Now, let us summarize the formula for de- It is also convenient to decom- 
pose de into the neutralino-exchange and chargino-exchange contributions as 
follows: 

These contributions have been extensively studied in the literature including 
in [26] which give 

4"' = -(j^EEMMJ^i^(3^) 

^ ' A=l X=l ^a: \ ex / 

where 

1 21nx\ 

A(x) = — 3 - X H 

^ ^ 2(1 V 1 -xy 

and 

, 1 / 2a;lnx\ 

B(x) = — {l + x-\ . 

^ ' 2(1 -x)2 V 1-x J 
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